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We report on the measurements of the superconducting order parameter in the nonmagnetic 
borocarbides LuNi2B2C and YNi2B2C. Andreev conductance spectra are obtained from nanoscale 
metalhc junctions on single crystal surfaces prepared along three major crystallographic orientations: 
[001], [110], and [100]. The gap values extracted by the single-gap Blonder-Tinkham-Klapwijk model 
follow the BCS predictions as a function of temperature and magnetic field and exhibit a small 
anisotropy. These observations are robust and reproducible among all the measurements on two 
different sets of LuNi2B2C crystals and one set of YNi2B2C crystals. We suggest possible explana- 
tions for the small gap anisotropy based on the recent detailed gap measurements by angle-resolved 
photoemission spectroscopy and the tunneling cone effect. Our results provide a consistent picture 
of the superconducting gap structure in these materials, addressing the controversy particularly in 
the reported results of point-contact Andreev reflection spectroscopy. 

PACS numbers: 74.20. Rp, 74.45. -l-c, 74.50.-l-r, 74.70.Dd 



I. INTRODUCTION 

Discovered more than fifteen years ago,^'^ the family of 
quaternary intermetallic compounds, RNi2B2C (R = rare 
earth elements), have attracted great attention<^"— They 
provide a unique opportunity for the investigation of the 
interplay between superconductivity and magnetism in a 
homologous series of compounds. They can be classified 
into four subgroups: magnetic members (R = Pr, Nd, 
Gd, Tb), superconductors (R = Sc, Y, Lu), coexisting 
materials (R = Dy, Ho, Er, Tm), and heavy fermions 
(R = Yb)/i^^— The non- magnetic superconducting boro- 
carbides, LuNi2B2G and YNi2B2C, a pristine model to 
study the superconductivity in this family, have the high- 
est superconducting transition temperature: Tc ~ 16.5 
K and 15.5 K, respectively. Despite intensive studies 
over the past decade, there is no consensus on their de- 
tailed gap structures. How magnetism and superconduc- 
tivity coexist or compete in the coexisting members re- 
mains a topic of considerable research interest since a 
systematic evolution of both phase transitions is observed 
over a comparable temperature range (1 — 20 K). On a 
broader context, the interplay between superconductivity 
and magnetism has become a central theme in the stud- 
ies of novel and unconventional superconductors includ- 
ing heavy fermions, cuprates, and the recently discov- 
ered iron-based superconductors. The magnetic and su- 
perconducting orders in RNi2B2G have different origins, 
arising from the localized 4f electrons on R'^"'" ions cou- 
pled via the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction;^ whereas the superconducting condensation 
is associated with the itinerant electron bands residing on 
the Ni2B2 layers.^ 

Complex Fermi surface (FS) structures with multi- 



ple sheets have been predicted from band structure 
calculations^'^ and confirmed by several experimental 
techniques fS^ii There are three bands crossing the Fermi 
level, all having electron-like characteriiii While the 
small ellipsoidal FS from the 19th band is known to con- 
tribute very little to the electronic density of states, both 
the 17th and 18th bands have been observed to comprise 
the major FS sheets although their topology is quite dif- 
ferent from each other. For the magnetic borocarbides, 
the 'cushion'-like FS from the 18th band, coming exclu- 
sively from the Ni 3d^2_y2 and Sd^y derived states, is 
not affected by the magnetic moments of the rare earth 
ions and, thus, the superconductivity originating from 
this band can survive under the development of magnetic 
order liiii^ 

The superconductivity in the non-magnetic borocar- 
bides RNi2B2C (R = Y, Lu) has been associated with 
both the 17th and 18th bands. Thus, a fundamen- 
tal question to be addressed is comprised of two parts: 
(i) whether the superconducting order parameter con- 
sists of multiple components in these bands; and (ii) 
what is the order parameter symmetry and anisotropy 
in momentum space. Early on, Shulga et alJ^ re- 
ported on the evidence for a multi-band nature of the su- 
perconductivity from the temperature dependence mea- 
surements of the upper critical field in YNi2B2G and 
LuNi2B2G. A variety of experimental investigations such 
as photoemission spectroscopy^*, Raman scatteringi^, 
thermal conductivity^, specific heat^^, and ultrasound 
attenuation^'^ point to a large anisotropy in the super- 
conducting gap functions. More detailed information on 
the gap anisotropy is reported from the measurements 
of thermal conductivity^^'^° and specific healj^ii^ as a 
function of magnetic field orientation. Glear four-fold 
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oscillations are observed in both kinds of measurements 
and interpreted as strong evidence for an anisotropic gap 
structure with gap minima located along the [100] and 
[010] directions. In particular, Izawa et al}^'^^ report 
the possible existence of point-like nodes from the obser- 
vation of diminishing oscillation amplitude as the field di- 
rection changes from the a6-plane to the c-axis. This may 
be consistent with an extremely large gap anisotropy: 
7gap = >100. It is noted that these measure- 

ments can te obscured by the anisotropy in properties 
other than the superconducting order parameter, making 
their interpretations non-trivial;^!^ This can be an issue 
particularly relevant to the borocabide superconductors 
since many properties are known to be anisotropic such 
as the FS^'"— upper critical field^ and vortex lattice 
structurcj^^"— 

Maki and coworkers^Si^ propose a hybrid s-fg wave or- 
der parameter, in which the g-wave component is added 
as an ansatz to simulate the point-node structure implied 
by experimental observations. Kontani^ point out that 
point-like nodes can be formed in a phonon-mediated s- 
wave superconducting order parameter if antiferromag- 
netic fluctuations exist on some parts of the FS con- 
nected by a nesting vector. Earlier, band structure cal- 
culations based on local density approximation predicted 
a nesting feature on the FS with nesting vectors Q ~ 
27r(0.5/a, 0) and 27r(0,0.5/a) in the basal planed This 
is directly observed from two-dimensional angular cor- 
relation measurements on LuNi2B2C using the electron- 
positron annihilation radiation technique^ Interestingly, 
an antiferromagnetic order in the magnetic borocabide 
compounds is observed to have an ordering wave vec- 
tor Q,„ « 27r(0.55/a,0,0)^"^ very close to the nest- 
ing vector. The phonon-mediated pairing has long been 
attributed to the superconductivity in the borocarbide 
family.—!^ For instance, Martinez- Samper et alM re- 
port on scanning tunneling spectroscopy measurements, 
in which features due to strong electron-phonon cou- 
pling reminiscent of a Pb tunnel junction are observed. 
Their extracted coupling constant is highly anisotropic 
due to soft phonon modes as observed in inelastic neu- 
tron scattering experiments 42r— According to the the- 
oretical arguments by Kontani)^ the nesting feature in 
the FS would generally weaken the electron-phonon cou- 
pling along Q. Thus, this scenario can explain the ex- 
perimentally observed anisotropic gap structure in the 
nonmagnetic borocarbide superconductors. 

Spectroscopy based on the measurement of Andreev 
reflection conductance has been adopted frequently as 
a measure to investigate superconducting order parame- 
ters (see refs. IH and H and references therein). Called 
point-contact Andreev reflection spectroscopy (PCARS), 
it is a simple and versatile technique, with which two 
groups have reported results on the nonmagnetic boro- 
carbide family. Disparate claims are made from measure- 
ments using virtually the same technique. Raychaud- 
huri et al^ observe a large anisotropy (7gap ~ 5) be- 
tween the gap values along the [100] and [001] direc- 



tions. They interpret this as evidence for the s-|-g or- 
der parameter with point nodes as suggested by Maki 
et al^ Later, they interpret this as evidence for multi- 
band superconductivity^^ Naidyuk and coworkers do 
not report such large anisotropy4i"— Between these 
two groups, the temperature dependences of the gap 
value along the [100] direction are strongly contradic- 
tory. Thus, our motivation for this work was, in part, 
to address this controversy with our own detailed gap 
structure measurement in the non-magnetic borocarbide 
superconductors.'^^ 



II. EXPERIMENTS 

The working principle behind point-contact spec- 
troscopy is that the energy dependence of quasiparticle 
scattering at an interface between two metallic electrodes 
is reflected as a non-linearity in the current-voltage char- 
acteristics. Analysis of such nonlinearity provides im- 
portant information on the scattering sources. When ap- 
plied to a superconductor, the most relevant scattering 
process is Andreev reflection, which is essentially a scat- 
tering off the pair potential. Thus, PCARS relies on the 
strong energy dependence of this scattering at a normal- 
metal/superconductor interface. In order to avoid com- 
plications due to other scattering processes (elastic or 
inelastic), a point-contact junction must fulfill require- 
ments for the ballisticity. Ideally, its size needs to be 
smaller than the electronic mean free paths of the elec- 
trodes. In practice, nanoscale metallic junctions can be 
made by bringing a sharpened metal tip into contact with 
a superconductor using fine mechanical adjustments. 

We perform systematic investigations of the gap struc- 
tures in the nonmagnetic borocabide superconductors 
YNi2B2C and LuNi2B2C. Differential conductance spec- 
tra along three major crystallographic directions are 
taken from two sets of LuNi2B2C single crystals of dif- 
ferent sources and one set of YNi2B2C single crystals. 
Crystals from different batches are used for each direc- 
tion in the set LuNi2B2C #1. For the sets LuNi2B2C 
^2 and YNi2B2C, crystals for all three directions are 
prepared from same batches. These single crystals are 
grown by flux method using Ni2B flux.— Bulk resis- 
tivity and magnetization measurements show that they 
are of high quality, as summarized in Table I. 

For a point-contact junction oriented along the c- 
axis, the as-grown surface of a single crystal is used, 
since its normal is along that direction. For in-plane 
junctions along [100] and [110] directions, crystals are 
embedded into low-temperature epoxy, cut, and pol- 
ished until atomic smoothness (rms roughness ~10 A) 
is achieved J^^i^ The desired crystallographic orienta- 
tions of the exposed surfaces are identified by X-ray 
diffraction, with 5 degree accuracy with the intended 
directions obtained. Since pristine crystals may con- 
tain a degraded surface layer where superconductivity 
is strongly suppressed, which is not favorable to PCARS 
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TABLE I: Summary of sample characteristics and PCARS results. RRR (residual resistance ratio) = _R(300K)/i?(rc,on), (Tc,on 
— onset transition temperature); rc,zcro: zero-resistance transition temperature; ATc = T^^on — Tc.zero; Aq: gap energy at 
zero-temperature; 7gap = Amax/Amin: maximum gap anisotropy; u = 1 — Aooi/Ano: weight for g-wave component (see text). 
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due to its surface-sensitive nature, the crystals are etched 
slightly in aqua-regia for about 20 seconds to expose fresh 
surfaces prior to making point-contact junctions. For the 
PCARS measurement, a crystal is cooled down to about 
2 K in the liquid helium cryostat, and then a sharp gold 
tip is moved to engage onto the crystal surface by a fine 
differential micrometer until a reasonable junction resis- 
tance is reached. 

The junction resistance, usually ranges from sev- 
eral ohms to tens of ohms. Applying the Wexler's 
formula'^^ using these Rj values along with known ma- 
terials parameters for borocarbides (po ^ 1-30 /iJl-cm 
and I ^ 500 A as reported in ref. flGh , we infer that the 
junctions are in the ballistic limit. What is more impor- 
tant to ensure the spectroscopic nature is that our mea- 
sured conductance spectra are reproducible as shown be- 
low and free from non-spectroscopic effects such as local 
heating. The differential conductance (G = dl/dV) as a 
function of voltage {V) is directly recorded by the stan- 
dard four-probe lock-in technique over wide temperature 
(1.6 — Tc) and magnetic field (0 — 9 T) ranges covering 
the whole phase space for the superconductivity in these 
materials. The local critical temperature in the junction 
area, , is determined by the temperature dependence of 
the zero-bias conductance. The measured is in good 
agreement with the bulk Tc determined by four-probe 
resistance measurement. This is an indication that the 
properties of the bulk rather than a possible degraded 
surface layer are measured in our PCARS experiments. 
For brevity, we simply refer to T-! as Tc in the following. 



III. SUPERCONDUCTING ENERGY GAP 

A. LuNiaBaC 

Figure [T] shows one set of G{V) curves for LuNi2B2C 
#1 crystals, taken as a function of temperature along the 
three major crystallographic orientations: [001], [110], 
and [100]. Note that at the lowest temperatures the An- 
dreev reflection amplitude and peak positions in the three 



directions are comparable. The data are analyzed us- 
ing the Blonder-Tinkham-Klapwijk (BTK) model^i^ as- 
suming a single isotropic s-wave gap. In this BTK model, 
there are three parameters: the superconducting gap (A) 
extracted from the BTK fit, the barrier strength (Z) be- 
tween the normal metal and superconductor, and the 
quasiparticle scattering parameter (F), which accounts 
for the smearing of a conductance curve due to short- 
ened quasiparticle lifetime or depairing effect caused by 
inelastic scattering or magnetic field. Here, we focus on 
the extracted superconducting gap and its temperature 
dependence. 

Figure [2] shows the temperature dependence of the su- 
perconducting gap. A, in the [001], [110], and [100] direc- 
tion, respectively. As shown, they all follow the standard 
BCS-like curve, yielding 2Ao/kBTc ~ 3.5, 3.8 and 3.6 in 
the weak-coupling limit. This is in good agreement with 
the results obtained from break junction experimentsi^ 
At the lowest measurement temperatures (^ 2 K), the 
[001] and [110] junctions show comparable superconduct- 
ing gap values with Aqoi = 2.4 meV and Ano — 2.6 
meV, respectively. The extracted gap value along the 
[100] direction is ~ 2.3 meV, giving a small anisotropy 
with the other directions (7gap ^1.13). As shown in Ta- 
ble I, the bulk Tc of this crystal is lower than those for 
the others. Its origin is not clear yet but the 2Ao/fcB7c 
value is consistent within the LuNi2B2C #1 set. Our 
gap values are similar to those reported by Bobrov et 
where the c-axis and a&-plane gap sizes from one- 
gap BTK fit are 2.25 and 2.55 meV, respectively. We 
note that their samples for a6-plane junctions are pre- 
pared without well-defined orientations. In our study, 
dozens of point-contact junctions have been measured 
and their conductance spectra are quite reproducible as 
demonstrated by the histograms in Fig. [3l which counts 
the occurrences for different gap values at the lowest mea- 
surement temperatures. Aioo exhibits a bit scattered dis- 
tribution although the most frequently observed value is 
2.3 meV. This behavior might be indicative of the tun- 
neling cone effect, which is most pronounced along the 
minimum gap direction, as discussed later. 

The superconducting energy gaps obtained from a 
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FIG. 1: (Color online) Temperature dependence of the differential conductance spectra for the point-contact junctions on 
LuNi2B2C #1 along (a) [001], (b) [110] and (c) [100] directions. These G{V) curves are obtained from junctions showing the 
most frequent gap values in Fig. 3. The data are normalized by the conductance at the negative maximum bias voltage in (a), 
(b), and (c). 
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FIG. 2: (Color online) Temperature dependence of the super- 
conducting gap in LuNi2B2C #1 extracted from an analysis 
using the single-gap BTK model along (a) [001], (b) [110] and 
(c) [100] directions. The symbols are experimental data and 
the solid lines represent the BCS prediction. 



single-band BTK analysis show slight deviations from the 
BCS curve in the low temperature region, as seen in Fig. 
2(b) and in other data in the following. We also note that 
in the low temperature region the best-fit curves are not 
as satisfactory as for a known isotropic single-band su- 
perconductor such as Nb. Bobrov et have reported 
similar observations in their PCARS on LuNi2B2C crys- 
tals and demonstrated that their G{V) curves can be 
better fit by assuming a continuous distribution of the 
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FIG. 3: (Color online) Histograms plotting the frequencies of 
the superconducting gap values measured from point-contact 
junctions on LuNi2B2C #1 along (a) [001], (b) [110], and (c) 
[100] direction. 



gap function with double maxima. They interpreted this 
as evidence for an anisotropic gap structure or multiband 
superconductivity in LuNi2B2C. Although such an anal- 
ysis is not an unreasonable approach considering recent 
developments in the field, we view it as marginal since it 
does not provide corroborative evidence for multiple or- 
der parameters. Note that strong evidence for multiband 
superconductivity in MgB2 does not come from improved 
conductance fitting but from clear features due to multi- 
ple gaps as observed in the original data.— 

The magnetic field dependence of the superconducting 
gap was studied at the lowest temperatures around 2 K, 
as shown in Fig. 01 The G{V) curves were taken with 



5 





[001] A A 


H(T) 


1 .25 









0.25 




T=2.21KjO iO* 






1.0 


1.00 




4.0 






7.0 


0.75 


(a) V/ 





1.50 
CD^1.25 
^ 1.00 

CD 

1.50 
1.25 
1.00 



[110] ft A 


H(T) 







T=1.72K O H 


0.3 




0.6 . 




2.0 




4.0 




7.0 


(b) y 





[100] . n 


H(T) 









0.5 


■T=2.34K ¥\ r\\ 


1.0 




2.0 




-5,0 




8.0 



-4 4 
Voltage (mV) 



0.0 



1.6 



■"o.s 



0.4 





■ 


[001] 




T 


[110] 






[100] 


(d) 








0123456789 
Magnetic Field (T) 



FIG. 4: (Color online) Magnetic field dependence of the 
differential conductance for the point-contact junctions on 
LuNizBzC #1 along (a) [001], (b) [110] and (c) [100] direc- 
tions. The data are normalized by the conductance at the 
negative maximum bias voltage in (a), (b), and (c). The 
gap A extracted from a single-gap BTK analysis is plotted 
in (d) and the quasiparticle broadening factor F in (e), as a 
function of the magnetic field. The solid lines in (d) and (e) 
represent best-fit curves using A = Ao(l — H/Hc2y^^ and 
F = Fo + co\/H, respectively. Here, Fo — 0.23 meV for all 
three junctions; co = 0.60 meV ■ T~^/^ for the [001] and [110] 
junctions, and Co = 0.45 meV • T"^''^ for the [100] junction. 



the field applied along the corresponding crystallographic 
direction. With increasing field, the conductance peak is 
suppressed and broadened with its position moving to- 
ward the zero bias. As shown in Fig|3] (d), the super- 
conducting gap decreases with increasing field, roughly 
following the BCS behavior, A = Ao(l - H/H^2)^/'^. No 
field-dependent signatures for multiple gaps are observed. 
The r parameter increases with \/H dependence, as plot- 
ted in Fig|l](e). 

Another set of LuNi2B2C crystals from a different 
source are prepared in the same way in order to check 
sample dependence. The GiV) curves and the extracted 
gap values are displayed in Fig. [5]as a function of temper- 
ature. Similar behaviors in the gap function to the previ- 
ous data set are observed, following the BCS prediction 
and showing a small anisotropy. We note that 2Ao/A:bTc 
values for these samples are systematically larger than 
those in the first set by ^ (3 — 9)% although their Tc 
values are comparable, as listed in Table. This means 
that they are in the intermediate coupling regime. What 
causes this difference between the two sets of LuNi2B2C 
crystals is not clear. 



2.5 
2.0 
1.5 
1.0 
0.5 

3.0 
^2.5 
« 2.0 

^1.5 

< 

Q.1.0 
(0 

O 0.5 







T(K) 






1.619 






3.53 
5.86 
8.30 










- 12.22 


16.16 



(a) 

[001] 



-4 4 
Voltage (mV) 
J 1 I L. 



1.6 

1.4 

S?^1.2 
> 

g-i.o 

0.8 



A A T (K) 




2.068 




3.17 




6.73 




^.11 








14.88 


^ 15.89 




4 6 8 10 12 14 16 
Temperature (K) 

FIG. 5: (Color online) Temperature dependence of the su- 
perconducting gap extracted from a single-gap BTK analysis 
for point-contact junctions on LuNi2B2C #2 along (a) [001], 
(b) [110] and (c) [100] directions in comparison with the stan- 
dard BCS curve. The insets plot the conductance data as a 
function of temperature. 



B. YNiaBaC 

In terms of electronic structure, there is not much dif- 
ference between YNi2B2C and LuNi2B2CZiii Although 
the maximum bulk Tc differs by 1 K, this can be ac- 
counted for by a small difference in the electron-phonon 
coupling constant or the density of states at the Fermi 
level. The superconducting gap has been reported to be 
anisotropic in both compounds. Nonetheless, we have 
carried out PCARS measurements on YNi2B2C single 
crystals in order to rule out the possibility that the gap 
anisotropy may be somehow enhanced in this compound 
than in LuNi2B2C. As shown in Table HI the RRR values 
of these crystals are comparable to those of LuNi2B2C 
crystals. 

The temperature evolution of the conductance spec- 
tra and the gap energy from point-contact junctions on 
YNi2B2C are shown in Fig. [B] and Fig. [71 respectively. 
Again, these data are analyzed with the single-gap BTK 
model to extract the gap values. They all follow the 
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FIG. 6: (Color online) Temperature dependence of the nor- 
malized difTerential conductance for the point-contact junc- 
tions on YNi2B2C along (a) [001], (b) [110] and (c) [100] di- 
rections. The data are normalized with respect to the negative 
maximum bias voltage. 



BCS curve although some deviations are observed in 
the low temperature region, similarly to those observed 
in LuNi2B2C. The averaged gap values obtained from 
dozens of contacts at the lowest measurement tempera- 
tures of ^ 2 K are given in Table ID The maximum gap 
anisotropy 7gap is 1.25, larger than in LuNi2B2C, but 
the gap is maximum along the [100] direction. This is 
in contrast to the case of LuNi2B2C, where the gap is 
maximum along the [110] direction. 

As in the case of LuNi2B2C, the extracted energy gap 
closes at the bulk resistive Tc in all three directions. This 
is a strong indication that our PCARS results represent 
bulk properties instead of a possible surface layer that 
is degraded during experiments. Theoretically, it is pos- 
sible that the small gap anisotropy may be caused by 
the averaging effect due to microscopic surface facets. In 
order to check for this possibility, we took conductance 
data on crystals with and without undergoing the chem- 
ical etching process after polishing. Similar results are 
obtained from both cases, ruling out such possibility. 
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FIG. 7: (Color online) Temperature dependence of the super- 
conducting gap in YNi2B2C along (a) [001], (b) [110] and (c) 
[100] directions in comparison with the standard BCS curve. 



IV. DISCUSSIONS 

As presented above and summarized in Table U we 
have observed a small gap anisotropy from both nonmag- 
netic superconductors LuNi2B2C and YNi2B2C. Maxi- 
mum gap anisotropy ranges between 1.12 and 1.25. This 
observation is reproducible and robust among different 
sets of samples and measurements. Our results are ap- 
parently contradictory with many recent reports in the 
literature claiming for a strong anisotropy in the gap 
structure, as discussed in Sec. Ijir— Below, we put for- 
ward possible explanations for this discrepancy, based on 
the Fermi surface topology and tunneling cone effect. We 
note that, by nature, thermal conductivity and specific 
heat in the superconducting state are sensitive to the 
gap minimum, whereas inherently anisotropic gap fea- 
tures could be smeared out in PCARS since the momen- 
tum of an injected electron can distribute over a wide 
range of angle with respect to the nominal surface nor- 
mal. We begin by reviewing the Fermi surface topol- 
ogy revealed by band structure calculations and recent 
angle-resolved photoemission spectroscopy (ARPES);^ 
and quantum oscillation experiments ^ii 

The claimed nesting structure on the 17th band 
FS has been observed by different experiments 
such as ARPES,— and electron-positron annihilation 
measurements^, where the nested part only occupies 
a small portion of the Fermi surface According to 
Kontani,'^^ antiferromagnetic fluctuations on the nested 
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parts of the FS could induce a point-node-like gap min- 
imum along the nesting wave vector. Recently, Baba 
et alM have reported on the detailed gap structures in 
YNi2B2C measured by high-resolution ARPES. Their re- 
sults show that the superconducting gap is observed on 
multiple bands (the 17th and 18th) with very different 
momentum dependences. More specifically, it is revealed 
that the gap is highly anisotropic on the 17th band with 
two minima, whereas the 18th band has a nearly constant 
gap. The two minima are 1.5 and 2.3 meV on different 
parts of the 17th band FS. Their momentum directions 
are identified as [100], consistent with other reportsji^Ti^ 
It was claimed that the points on the FS where the gap 
shows a minimum of 1.5 meV can be connected by the 
nesting vector Q ^ 27r/a(0.55, 0, 0), implying an intimate 
connection between the minimal gap and the FS nesting. 
It was also argued that the reason why the gap minimum 
has a non-zero value is because their ARPES data were 
taken off the basal plane with ~ 0.5^. 

Based on these observations from ARPES, it is natu- 
rally understood that the node-like gap minimum along 
[100] cannot be observed in PCARS since the shape of 
a conductance curve will be dominated by the other gap 
minimum (2.3 meV) which is located along the same di- 
rection. Indeed, this gap value is close to those we ob- 
served in this study (see Table HI), explaining why we 
observe only a small gap anisotropy in both LuNi2B2C 
and YNi2B2C. 

Our experimental results are in strong contradiction 
with those reported by Raychaudhuri and coworkers^^*^ 
From PCARS measurements on YNi2B2C, they reported 
a gap anisotropy of 7gap ~ 4.5, much larger than ours 
(1.12— 1.25). This discrepancy is mainly associated with 
their extremely small gap value along [100] (~ 0.415 
meV) since their gap value along [001] (~ 1.8 meV) is not 
much different from ours. Moreover, they observed that 
Aioo closes at '-^ (7 — 9) K, much lower than the bulk Tc- 
This is in strong contrast with our observation of Aioo 
closing at the bulk Tc. We note that Naidyuk et al.'^'' also 
reported similar temperature dependences to ours. Ray- 
chaudhuri et al. originally interpreted their results as evi- 
dence for s-l-g gap structure with point nodes^^ and later 
on claimed for the multiband nature of the supercon- 
ductivity in borocarbides.'*^ We think this is an unlikely 
interpretation considering the recent ARPES measure- 
ments discussed above since there is no way for PCARS, 
which doesn't have such high momentum resolution as 
in ARPES, to probe only an extremely small gap value 
0.415 meV) while there exists another but sizable gap 
minimum (^ 2.3 meV) along [100]. Rather, it is more 
natural to attribute the exotic behavior of their measured 
Aioo to non-intrinsic effects such as suppressed supercon- 
ductivity on the sample surface, as was also pointed out 
by Naidyuk et al^ 

Next, we discuss the tunneling cone effect on the con- 
ductance spectra. In the original BTK model,- - the in- 
tegration is carried out only over the energy, assuming 
an isotropic gap structure. Thus, it is inaccurate for the 



analysis if the energy gap has substantial dependence on 
momentum. A straightforward way to deal with this sit- 
uation is to include an integration with respect to mo- 
mentum. Here sets in the tunneling cone effect. '''^'^^ If we 
consider a planar tunnel junction with a potential barrier 
of finite thickness, the transmission probability varies de- 
pending on the momentum direction, being maximal for 
a direction normal to the barrier. Assuming a Gaus- 
sian distribution for the probability, the tunneling cone 
angle can be as small as 5 — 10 degrees in typical tun- 
nel junctions)^ This angle is expected to increase with 
decreasing barrier strength. Therefore, its effects are 
not negligible in PCARS due to inherently small barrier 
strength as well as small junction size (uncertainty prin- 
ciple). Here, we simulate conductance spectra in order 
to study the effect of tunneling cone on the gap structure 
with large anisotropy. Namely, we calculate conductance 
curves that are obtained from a given gap function with 
a varying cone angle, and then determine the gap values 
from the peak positions. 

In order to simulate conductance curves, we adopt 
the extended BTK model for a d-wave superconductor, 
formulated by Tanaka and Kashiwayaji^ but only con- 
sider gap functions without a sign change since there is 
substantial evidence in these materials for this premise. 
Then, the conductance kernel <ts{E) for a given ^^r, an- 
gle of incidence in the normal-metal electrode, can be 
written as: 



as{E) = (TAT 



i + ajv|r+|^ + (aAr-i)|r+r_ 
|i + ((Tjv-i)r+r_|2 



(1) 



where cjn = Z = and r± = ^ . 

Here Zq is the barrier strength, and A± represent the 
pairing potentials for the transmitted electron-like quasi- 
particles or hole-like quasiparticles, respectively. Thus, 
for a normal-metal/superconductor junction, the total 
conductance axiE) is given by the integration of as{E) 
with respect to a solid angle ft, 



<Jt{E) 



J dVLas{E) cos 6'jvP(6i at) 
j dilaN cos 9 iqP [On) 



(2) 



where P{9n) is the transmission probability for a given 
On- As an approximation, we assume a Gaussian-type 
distribution: 



P{eN) oc e 



(3) 



The parameter 8d is a variable characterizing the tun- 
neling cone for a given junction. We consider two 
gap functions: s-l-g wave as proposed by Maki et alM 
and anisotropic s-wave gap, both having four-fold point- 
nodes. The s-|-g gap function can be written as A(k) — 
iAo[l — sin'' 9cos{4:4>)] and for an anisotropic gap we as- 
sume A(k) = Ao[l — sin'' cos2(20)] for simplicity. In 
the zero-temperature limit, the peak position of a con- 
ductance curve can be taken as the energy gap. Fig- 
ure [5] plots thus obtained A along three major directions 
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FIG. 8: (Color online) Dependence of the peak position, 
Apeafe/Ao, of the conductance curve simulated by the ex- 
tended BTK model on the tunneling cone angle, Od, assum- 
ing a nodal gap structure of (a) s-|-g wave and (b) anisotropic 
s-wave, respectively. 



as a function of the tunneling cone angle, Qd- An as- 
sumption is made that the tunneling cone Qd does not 
depend on the surface orientations. The inherently large 
anisotropy is observed for a narrow tunneling cone. How- 
ever, the anisotropy becomes smaller as the cone angle 
increases. Aioo = Aqoi at S/j —25 and 35 deg., for 
the s-|-g and anisotropic s-wave case, respectively. Al- 
though qualitatively similar, there is difference between 
the two gap functions in their response to the tunnel- 
ing cone. That is, the calculated gap values do not 
merge in the s-|-g case, whereas they merge into one 
for large tunneling cone in the anisotropic s-wave case. 
This indicates that the anisotropic s-wave with Qd > 30 
deg., rather than the s-|-g, agrees with the small gap 
anisotropy we have observed constantly. If we consider 
more general case of the s-|-g gap function, written as 
A = Ao[(l — w) — sin^ 6* cos(4(/))], we can estimate the 
weight for the g-wave component. It turns out that uj is 
very small in all our measurements as listed in Table ID 
Thus, we argue that the s+g wave is unlikely as a candi- 
date gap function in borocarbides, in contrast to several 
reports in the literatureJ^i^Sii^ 

Our simulation based on tunneling cone effect is essen- 
tially similar to the BTK fitting with a distribution of 
gap values as reported by Bobrov et alM- in that contri- 
butions from different gap values are added in both cases. 
As we do not claim for an anisotropic gap structure in 
the borocarbide superconductors from our tunneling cone 
simulation, we don't view their arguments as corrobora- 
tive evidence for the multiband superconductivity, either. 

Studies on the thermal conductivitj*^ and specific 
heat*'^ in doped single crystals of Y(Nii_a;Pt2:)2B2C have 
shown that the gap anisotropy is washed out as the 
doping level increases. This was interpreted as due to 
gap opening induced by the nonmagnetic impurity, in 
agreement with theoretical calculations. ^^'^^ Our trans- 
port measurements indicate that the samples used in our 
PCARS are of high quality and no correlation between 



RRR and ^gap is observed, as summarized in Table ID We 
note that the Aioo value of a LuNi2B2C #2 single crys- 
tal changes very little after annealing although its RRR 
increases to 34.3 (not shown here). Moreover, the spe- 
cific heat results reported in refs. 21i and i22i which were 
obtained using single crystals from the same sources as 
ours, have shown clear four-fold oscillations under mag- 
netic field. Thus, the possibility of impurity scattering 
being an intrinsic origin for the small gap anisotropy can 
be excluded in our PCARS study. 

Considering our experimental observations in relation 
to the discussions presented above, we predict that the 
anisotropic gap structure can be better probed by thin- 
film based tunneling spectroscopy. In particular, the two 
gap minima along the [100] direction as seen in the re- 
cent APRES study^^ should also be discernible in the 
tunneling conductance data. 



V. CONCLUSIONS 

We have carried out PCARS measurements on the non- 
magnetic borocarbide superconductors LuNi2B2C and 
YNi2B2C in order to investigate their detailed gap struc- 
tures. Conductance spectra as a function of crystallo- 
graphic orientation, temperature, and magnetic field are 
taken from two different sets of LuNi2B2C crystals and 
one set of YNi2B2C crystals. Analysis based on the 
single-gap BTK model shows that the superconducting 
gap follows the BCS predictions as a function of tem- 
perature and magnetic field, closing at bulk supercon- 
ducting transitions, albeit a slight deviation is sometimes 
observed. The measured superconducting gaps exhibit a 
small anisotropy among the three major directions: [001], 
[110], and [100]. Any clear evidence for multiple gaps is 
not observed. These results are completely reproducible 
and robust in all our measurements. In order to explain 
our observations, two possible scenarios are considered in 
terms of the recent ARPES studies and tunneling cone 
effect. We argue that these two effects may render it im- 
possible to observe the inherently large gap anisotropy 
using PCARS. Thin-film tunneling spectroscopy on the 
[100] surface should better detect the gap minima, which 
then can serve as a confirmation of the anisotropic struc- 
ture and multiband nature of the superconducting or- 
der parameter in these compounds. Our study addresses 
the controversy over the gap structure in the nonmag- 
netic borocarbides, particularly among the studies using 
PCARS techniques. 
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